(N ' b 



UCLA/00/TEP/27 

Alternative Dimensional Reduction via the 

Density Matrix 

C. A. A. de Carvalho a ' f f] and J. M. Cornwallf] 

q | a Department of Physics and Astronomy, University of California, Los 

Angeles, CA 90095-1547 
Permanent address: Instituto de Fisica, UFRJ, CP. 68528, Rio de 
Janeiro, RJ 21945-970, Brasil 
A. J. da Silvafl 
Instituto de Fisica, Universidade de Sao Paulo, 
CP. 66318, Sao Paulo, SP 05315-970, Brasil 

Abstract 

We give graphical rules, based on earlier work for the functional 
Schrodinger equation, for constructing the density matrix for scalar 
and gauge fields in equilibrium at finite temperature T. More useful 
is a dimensionally-reduced effective action (DREA) constructed from 
the density matrix by further functional integration over the argu- 
ments of the density matrix coupled to a source. The DREA is an 
O 1 effective action in one less dimension which may be computed order 

by order in perturbation theory or by dresseddoop expansions; it en- 
codes all thermal matrix elements. We term the DREA procedure 
alternative dimensional reduction, to distinguish it from the conven- 
^ ■ tional dimensionally-reduced field theory (DRFT) which applies at 

infinite T. The DREA is useful because it gives a dimensionally- 
reduced theory useable at any T including infinity, where it yields 
the DRFT, and because it does not and cannot have certain spurious 
infinities which sometimes occur in the density matrix itself or the 
conventional DRFT; these come from InT factors at infinite temper- 
ature. The DREA can be constructed to all orders (in principle) and 
the only regularizations needed are those which control the ultraviolet 
behavior of the zero-T theory. An example of spurious divergences 
in the DRFT occurs in d=2+l 99 4 -theory dimensionally reduced to 
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d=2. We study this theory and show that the rules for the DREA re- 
place these "wrong" divergences in physical parameters by calculable 
powers of InT; we also compute the phase transition temperature of 
this </? 4 -theory in one-loop order. Our density- matrix construction is 
equivalent to a construction of the Landau-Ginzburg " coarse-grained 
free energy" from a microscopic Hamiltonian. 

1 Introduction 

The study of relativistic field theories at finite T and/or density 0, in 
or out of equilibrium ||, is quite mature but continues to be of great inter- 
est, for example, in connection with RHIC and CERN relativistic heavy-ion 
collision experiments searching for the QCD transition from a hadronic to 
a quark-gluon plasma phase; a possible color-flavor locking transition to a 
superconducting phase of QCD at finite particle density; the transition to a 
disoriented chiral condensate of pions that could provide an explanation for 
the so-called Centauro events observed in cosmic rays; and cosmological tran- 
sitions [f| . There is an equally-extensive body of literature on non- relativistic 
systems at finite T and so forth, the preponderance of which deals with crit- 
ical phenomena ||. 

In this paper we will be concerned with the description of all these systems 
via the density matrix and a related dimensionally-reduced effective action; 
the methods used can be easily extended to the closely-related Wigner distri- 
bution functionfj We begin by finding rules for describing the density matrix 
as the negative exponent of a dimensionally-reduced action Sp. These rules 
follow from a microscopic theory by the usual procedure of integrating over 
Euclidean time evolution. This action, however, is not our final goal, because 
it may contain spurious ultraviolet divergences. If spurious divergences do 
occur, they will also occur in the usual dimensionally-reduced field theory 
(DRFT) arrived at by the standard procedure |J of going to infinite T. So 
while it is interesting to know how to construct the density matrix system- 
atically, it is of less use than the ultimate result of our procedure, which is a 
dimensionally-reduced effective action (DREA) constructed by further func- 

The density matrix is a diagonal matrix element; the Wigner distribution function is 
effectively an off-diagonal matrix element. The reader will see that it is quite straightfor- 
ward to extend the density-matrix rules to the Wigner distribution function; we will not 
discuss the Wigner distribution function rules in detail here. 
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tional integration over the density matrix, coupled to sources. The DREA 
has no spurious divergences and contains in it all possible thermal matrix 
elements which could be constructed by conventional means. The procedure 
adopted here is a natural extension of a construction used JF| in the functional 
Schrodinger equation. 

In the case of d=0+l theories (quantum mechanics) in many situations 
the integral over the Euclidean time evolution can be done explicitly f|, ||, [K| , 
yielding integrals which are complicated functions of T; in field theory it 
is, of course, impossible to do the corresponding Euclidean time integrals 
exactly. But one must expect similarly-complicated dependences on T in 
the field-theoretic effective action, and these dependences will be important 
in estimating such properties of thermal field theories as phase transition 
temperatures. 

For relativistic systems it is common to use the DRFT, which is appro- 
priate in principle only for infinite T, at temperatures very high compared to 
typical energies of the zero-T field theory. The DRFT may (we will study an 
example) have divergences which are not present at any finite T; while this 
does not impair the essential correctness of the dimensionally-reduced the- 
ory it does introduce arbitrary parameters which should not be there. The 
point is that at any T, however large, a field theory may have no ultraviolet 
divergences not already contained in the zero-T field theory; this well-known 
consideration follows from the Feynman rules for finite-T field theory. But 
sometimes a dimensionally-reduced field theory does have "wrong" diver- 
gences (which must, in fact, be absent in physical quantities) which interfere 
with the use of the DRFT as a good approximation for large but finite T. 
An example is finite-T d—2+1 ip A theory, where the dimensionally-reduced 
theory is just d=2 (p A theory. This latter theory has a one-loop logarithmic 
mass divergence which is absent in the d=2+l theory, and we will give rules 
for the density matrix for the d=2+l theory which automatically remove this 
divergence (replacing it by a finite and calculable factor of In T) . 

For nonrelativistic systems there are corresponding problems of renormal- 
ization, but they are generally bypassed in the study of second-order phase 
transitions. In studying such transitions, the need for ultraviolet renormal- 
ization and the fact that one has to start from a microscopic theory can 
be relegated to secondary status, because one is concerned only with long- 
range properties. Most treatments of nonrelativistic systems exhibiting phase 
transitions, such as in condensed matter physics, are not concerned with a 
detailed microscopic knowledge of the Hamiltonian, but rather with a uni- 
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versal description via the renormalization group that captures the essence 
of the long-range collective behavior responsible for the transition and so 
renormalization is not intrinsically required. The use of the renormaliza- 
tion group in such systems is a tool to extract the long-range behavior: it 
provides a way of eliminating short-range scales in favor of collective long- 
range parameters by integrating out short- distance scales to obtain effective 
Hamiltonians, and investigating the appearance of fixed points associated 
with criticality. Near those fixed points, it allows for the calculation of the 
exponents that characterize the critical behavior of physical quantities. The 
theories studied can be (and often should be) defined on a lattice; whether 
this is the correct microscopic description or not makes no difference to the 
long-range properties. In other words, the tools used to study long-range 
behavior are not capable of calculating parameters such as phase-transition 
temperatures, whose values depend on phenomena at shorter distances and 
in particular on the detailed functional dependence of the Hamiltonian on T. 
This Hamiltonian is typically a Landau-Ginzburg phenomenological "coarse- 
grained free energy", whose parameters' dependence on T and, possibly, on 
external parameters is chosen purely phenomenologically; it is not derived 
from a fundamental picture. 

This article addresses the question of how actually to derive the DREA, 
starting from the microscopic Hamiltonian that describes a field theory. This 
question is essentially equivalent to asking how to find the Landau-Ginzburg 
coarse-grained free energy, as the dimensionally-reduced theory whose ma- 
trix elements are those of the DREA itself. One can thereby replace the 
phenomenology of the Landau-Ginzburg free energy by explicit calculation 
of, say, the T-dependence of this free energy. In the case we study in this 
paper the underlying microscopic theory is considered to be known (a rel- 
ativistic quantum field theory), and the density matrix and the DREA are 
to be found. We will consider scalar and gauge field theories; nothing much 
new is found by studying fermionsf]. Just as with any effective action, the 
DREA contains in its infinitely-many terms all the (one-particle irreducible) 
matrix elements which can be constructed by tracing time-independent oper- 
ator products with the density matrix. But actually constructing the density 
matrix itself is not as useful, since this object contains the same spurious 
infinities that the DFRT contains. 

5 Except possibly for topological effects involving fermions, which we will not discuss 
here. 
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The technical reason why these spurious infinities occur is that the prop- 
agators to be used in integrating over Euclidean time dynamics to form the 
density matrix are not the standard thermal propagators, and these non- 
standard propagators do lead to spurious ultraviolet divergences. However, 
in constructing the partition function from the density matrix, the density- 
matrix propagators combine with integration over certain source terms in- 
volving the time-independent fields which are the argument of the density 
matrix and cancellations occur; the result is equivalent to using the usual 
thermal propagators. As a result, the partition function can only involve 
ultraviolet divergences (and subsequent renormalization) of the T = d + 1- 
dimensional theory. The same combination of propagators and source terms 
occurs in constructing the DREA, which therefore has only the T = diver- 
gences. 

We note here some general properties of the DREA: At T large compared 
to any mass scale m in the original zero-temperature theory, local terms[] 
in the DREA with large powers of its argument 0(x) will be accompanied 
by large negative powers of T and vanish at infinite temperature; this is 
the way that local terms in the DREA yield the DRFT. As T — * 0, there 
will also be factors like exp(— m/T) along with these powers of T, and the 
zero-temperature limit is well-defined. In this limit the DREA is precisely 
the (negative logarithm of) the ground-state Schrodinger wave functional 
0. The point which is of particular interest in the present paper is the 
appearance of marginal terms depending on In T which lead to the unwanted 
divergences in the DRFT at infinite temperature. 

The paper is organized as follows: Section ^| defines the density matrix 
for a scalar field theory, shows how to compute the reduced theory perturba- 
tively, discusses how its ultraviolet behavior is related to the renormalization 
of the original theory, and outlines the construction of the DREA. Section ^ 
discusses some general questions of gauge invariance of the DREA for gauge 
theories. Section ^ illustrates how the dimensional reduction obtained via 
density matrices can be used in practical applications: a modified reduced 
theory (defined by its Feynman rules) is proposed, which should yield the 
same results as the one derived from the density matrix, and a perturbative 
discussion of the phase transition for a scalar theory in two spatial dimen- 

6 Local terms in an effective action are candidates for the action of a field theory which 
could have the DREA as its effective action. Non-local terms correspond to one-particle 
irreducible Green's functions. 
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sions is presented. This last application is intended to emphasize that our 
approach reduces the discussion of the phase transition to the study of an 
(almost) ordinary system in d dimensions, plus the d+1 subtractions required 
by the usual renormalization procedure. 

2 The density matrix for scalar theories 

The partition function for a self-interacting scalar field theory in contact with 
a thermal reservoir at T (f3 — 1/T) can be written as a functional integral 
over the density matrix p: 

Z(0) = jpD0]p[/3;0,0], (1) 

p[/3;0,0] = ^]e- 5M , (2) 
S[<p] = dr J rf d x (£ F [<p] + d[<p]) , (3) 

C F [ V } = \{(d TV f + iy V f + m 2 V 2 } . (4) 

Here Hp and £/ are, respectively, the free and interacting Lagrangeans. We 
will be interested in <p 4 interactions (both broken and unbroken). The field 
0(x) is the boundary value of ip(r, x) at both r = and r = j3, that is, the 
integral § is to be performed over all ip that satisfy the boundary conditions 
ip(0, x) = <p(/3,x) = 0(x). The density matrix is a functional of <fi only. 
The remaining integral over the 0's is unrestricted (except for the vacuum 
boundary conditions that must be imposed at spatial infinity). 

We may write p = e~ s °, Sp being a certain T-dependent dimensionally- 
reduced action. This action is NOT the DREA, whose construction needs 
further discussion. The field which is the argument of Sd depends only on 
the d spatial coordinates; all the r dependence of the original d + 1-theory 
has been eliminated through the tp integration. The fields 0(x) are the nat- 
ural degrees of freedom of the reduced theory. Any thermal observable can 
be constructed by integrating over the fields 0(x) an appropriate function 
of <fi weighted with the corresponding diagonal element of the density ma- 
trix. The Euclidean time evolution can be viewed as an intermediate step 
which calculates the weights. In carrying out the functional integrals over 0, 
one would notice certain cancellations between divergences occurring in the 
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construction of Sp and divergences in the functional integral. These cancel- 
lations remove the spurious ultraviolet divergences (those not occurring in 
the T = theory). 

To complete the construction of the DREA it is only necessary to in- 
troduce further sources J(x) coupled to the field 0(x) through the usual 
term / J<p. The logarithm of the resulting functional integral over <p is then 
Legendre-transformed to yield the DREA. 

Note that neither the density matrix nor the DREA is, in general, of 
the form e~^ H with H being independent of f3 as in ordinary statistical 
mechanics; their (3 dependence is far more complicated (which justifies the 
"almost ordinary" of the previous Section). This had already been pointed 
out in the analogous discussion of the transfer matrix carried out in Ref. 
[TT| . The density matrix provides a direct but alternative way of deriving a 



dimensionally-reduced theory. We will proceed to construct it perturbatively, 
and indicate where appropriate the generally straightforward generalization 
to dressed-loop expansions. 

2.1 The perturbative construction for scalars 

The integral over the fields (p in equation (H) has to respect the boundary 
conditions. For the free theory, the integral is quadratic and can be computed 
exactly by the saddle-point method. We must solve the free field equation of 
motion, subject to the boundary conditions: 

(-<9 T 2 -V 2 + m 2 Mr,x)=0, (5) 

^(O,x) = ^O3,x)=0(x). (6) 

Fourier transforming in x leads to an ordinary differential equation. The 
solution satisfying the boundary conditions is: 

rzu ^\ f j3l c °sh[w k (r - (3/2)} ~ - k . x 

^ x) = J cosh ( /W2) m e (7) 

f '"' k 0(k) e* k ' x , (8) 



(27r) d 

where = +[k 2 +m 2 } 1 ^ 2 . It depends functionally on 0(x). Its r dependence, 
however, is completely specified. 
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We now expand tp around (p, treating it as a thermal background for the 
free theory: 

<p(r, x) = tp(r, x) + r/(r, x) , (9) 

7j(0,x) = »K/3,x) = 0. (10) 
(Similarly, the <p propagator is the sum of the (p propagator and the 77 propa- 
gator.) As usual, the fluctuation 77 has to vanish at r = and r = f3 because 
(p already satifies the boundary conditions. The expression for the density 
matrix element becomes: 

p\j3; 0, 0] = e~ SD = e- SF ® £[Dr]}e~ SF[v] - Sl & +v] , (11) 

where Sf and Si refer to free and interacting actions, respectively, and the 
path integral sums over fluctuations that vanish at r = and r = f3. 

Note that equation ([11]) is exactly of the same functional form as given for 
the functional Schrodinger equation in Ref. 0; it is a d + 1-dimensional func- 
tional integral (with an extra factor exp —Sp[0\), with sources depending on 
(p through the interaction term Si[(p + 77]. The free propagators used to eval- 
uate the functional integral differ from those of the functional Schrodinger 
equation, but the structure is exactly the same, and so the remarks made in 
concerning the extension of perturbation theory to dressed-loop expan- 
sions apply without change to the present case. We therefore need not repeat 
those remarks here. 

For the free theory, Si = and we can compute the integral exactly 
to obtain the fluctuation determinant (to power —1/2) using the techniques 
of Refs. m H [UJ. For the interacting theory, we first perform a functional 
Taylor expansion of the interacting action around the thermal background, so 
as to obtain a series (or a polynomial in the ip 4 case) in powers of 77. Then, we 
expand the exponential e~ Sl in a power series. The integral over fluctuations 
will Wick-contract the various products of 77's which appear multiplied by 
derivatives of the interacting Lagrangean taken at the thermal background. 
As a consequence, one is led to compute: 

(r7(r 1 ,x 1 )---77(r fc ,x fc )) = A F 1/2 ^[I?77]e-^ M 77(r 1 ,x 1 )---7 / (r fc ,x fc ) . (12) 

Af stands for the fluctuation determinant of the free theory, a (infinite) 
normalization constant that depends on f3. The result is simply: 

(77(7-1, xi) • ■■v{r k ,x k )) =J2 G v( T h, x h;n 2 ,Xi 2 ) ■ • •G r) (r ifc _ 1 ,x ife _ 1 ;r ife ,x Jfe ) , 

(13) 
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if k is even, and zero otherwise. J2p denotes sum over all possible pairings of 
the {ij}. There remains to compute the propagator G v , which must satisfy: 



'-d 2 T - V 2 + m 2 )G v (r, x; r', x') = 8{r - r')<5 d (x - x') 
G„(0,x; r / ,x') = G r? (Ax;r',x')=0, 



(14) 
(15) 



because of the vanishing boundary conditions for the fluctuations. If we 
Fourier transform the spatial coordinates, we obtain an ordinary differential 
equation which leads to J/J |9], |10[ : 



GJr, x; r',x' 



d d k sinh(w k r<) sinh[w k (/3 - r>)] e i k .(x-x') 
(27r) d U7 k sinh(/5wk) 



(16) 



where T < (r > ) = min(max){r, r'}. Note that this propagator is not r-translation 
invariant. Formally, we have: 

oo 

p\P\ <f>, 0] = Polfl; 0> 0] + E Pm[/5; 0] , 



Po 



m=l 



Pr, 



PO 



ml 



E 



(n 

=0 j=l 3 



-1/2 



n 



(r] ni (l)---T] nm (m)) 



j ■ 



(17) 

(18) 
(19) 



where Jl- = c?Tj J c/ d x.,-, the argument of 77 (j) stands for (r 3 -, x 3 -), and £j 
is the rij-th derivative of £/, computed at the thermal background (p. Ex- 
pression (|T7j) depends functionally on through </?. 

The spatial integrations should be converted to momentum space, where 
the propagator is diagonal (oc 5 d (\t' — k)). Performing the integral over k will 
lead to ultraviolet divergences. Therefore, we introduce a momentum cutoff 
A. It will be convenient to think of such a cutoff as analogous to an inverse 
lattice spacing in a lattice regularized version of the theory. In the cutoff 
theory, we can compute each term in the series and obtain the T and cutoff 
dependence of the various coefficients of a functional expansion in 0(x) (or 
0(k)). The logarithm of ( |T7D will then be —Sd, a cutoff action which will 
enter further functional integrals over 0(x). 



9 



2.2 The ultraviolet behavior 



The construction of the previous section allows us to find a cutoff reduced 
action Sd order by order in perturbation theory. We will now specialize our 
discussion to Cj = Ay? 4 /4!. As A — ► oo, the reduced theory exhibits the 
ultraviolet divergences of both d + 1 and d dimensions (the Appendix shows 
this explicitly for d = 2). When computing correlation functions, the former 
are to be eliminated through the usual renormalization procedure, which 
requires the same subtractions as in the T = 0, d + 1-dimensional, original 
theory. The d- dimensional divergences, however, are a consequence of failing 
to complete at the density-matrix stage the two-step process of integrating 
over the fields T](t, x) followed by integration over the 0(x) fields which are 
the arguments of the density matrix. The ^-generated divergences should 
automatically cancel upon doing the final integration over the remaining 
0(x) which is necessary to find the partition function Z(f3). Indeed, this is 
what happens, as we will show. There is no reason that these divergences 
should be absent from Sd itself. 

We can find these cancellations by calculating Z(j3). After having done 
the integration over the fluctuation fields rj, we still have to integrate (in 
momentum space) over the remaining 0(p) in order to compute the vacuum 
bubbles in the expansion ([H]). The free part of the reduced action Sd for 
0(p) is given by Sp[0\- An integration by parts, and use of the equation of 
motion satisfied by (p lead to: 

SfW\ = \J ^0(p)[2 Wp tanh(/3^ p /2)]0(-p) , (20) 

which defines the free 0-propagator as: 

G*(p,p') = (27r) d 6 d (p + p')[2w p t a nh((3w p /2)]- 1 . (21) 

This propagator will be used to contract any two <fi fields. In the language 
of Feynman graphs, any line representing a contraction of the corresponding 
fields appearing in (|TTD will have a factor: 



coshMr- ^/2)] [2u;ktanh( ^ k/2)rl ooshMr' - (3/2)} 

cosh{f3w k /2) [ KH 1 n cosh(/3w k /2) v ; 

For contractions of r) fields, we must use the G v propagator. Any line repre- 
senting such a contraction will have a factor: 

sinh(wkT') sinh[u> k (r — /3/2)]©(r — r' 



w k sinh(/?w k ) 



+ [t <-> t] . (23) 
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Note that this not the usual thermal propagator. However, the combination 
of the integration over rj and together leads to a cancellation of the ultra- 
violet divergences appearing in the separate integrals, as we see by rewriting 
the products in terms of sums of hyperbolic cosines of arguments (V — r) 
and (V + r + (3). The dependence on the latter cancels when we add (|2~2"D to 



23|). The sum may be written as [12 



^ , , _ k] = cosh{w k [/?/2 - |t'-t|]} If* cos[27rnT(r / - r)] 
thlT T ' J " 2 Wk sinh(/W2) /3 n f^ (27rnT)2 + w 2 ■ 

(24) 

It is worth noting that the first form of G t h can be written in a familiar form 



which Pisarski[[y|] has advocated for its practicality in calculations: 



„ 1 P — ivu\t'-t\ p Wk\r'-T\ 

GtJr'-r-k) = [ 5 2—1. (25) 

v ' 7 2w k L l-e-^k i_ e /W v ; 

Clearly, adding the rj and lines reconstructs the thermal (finite T) propa- 
gator of the d + 1-dimensional theory, which is r-translation invariant. Every 
graph in the expansion in vacuum bubbles of the original theory with the 
propagator (pi]) may have each of its lines replaced by the sum of a //-line 
plus a 0-line. This will generate all possible combinations of rj and 0-lines 
which appear the (cutoff) reduced theory. This proves that integrating the 
reduced theory over cf) will reproduce the results of perturbation theory to 
all orders, and that the only subtractions that are required are those of the 
T = theory in d + 1 dimensions. 

The final step to the DREA is very easy: Instead of constructing the 
unadorned partition function by functional integration, as in equation ([I]), 
we construct the partition function in the presence of a current J(x): 

Z[f3- J] = exp(-W[(3; J}) = J V<f>exp[S D + J d d xJ(x)(j)(x)}. (26) 

Now all the integrations over both rj and <fi have been done, so the same 
cancellations as in the partition function itself still occur. One then only 
needs to Legendre-transform from W[j3; J] to the effective action r[/3; (0)]; 
this, the DREA, is free of spurious ultraviolet divergences. 

Before applying these considerations to a specific field theory, we discuss 
how the general principles apply to gauge theories. 
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3 Density matrix for gauge theories 



Begin by introducing notation. We use the standard anti-Hermitean matrix 
form of the gauge potentials, with the coupling constant g absorbed in the 
potential: 



Here the A a /2 are the group generators in conventional normalization. Gauge 
transformations act on via the operation of unitary matrices U: 



It will sometimes be convenient to denote the time component of the d + 
1-dimensional potential Ao as an adjoint scalar (in d dimensions) field 
Greek indices refer to the d + 1-dimensional theory and Latin indices to the 
dimensionally-reduced theory. The potential A^ depends on both r and x. 
There is a corresponding potential depending only on x which is the 
argument of the density matrix or of the DREA; this can also be written as 



The extension of our formalism to gauge theories requires some attention 
to questions of gauge invariance and gauge fixing. The partition function 
can be written as the trace of exp(—/3H) where H is the microscopic time- 
independent Hamiltonian of the gauge theory, and the set of basis functions 
for the trace can be taken as the eigenfunctions of the functional Schrodinger 
equation. Since H commutes with the generator of time-independent gauge 
transformations, this basis may be chosen to be invariant under such gauge 
transformations.^ Therefore all matrix elements (of gauge-invariant operator 
products) found by tracing over the density matrix are gauge-invariant. But 
this does not necessarily mean that it is convenient to calculate or present 
the density matrix as completely gauge-invariant; it proves more convenient 
to present it in a form where there is a gauge-fixing term. This c/-dimensional 
gauge-fixing term is inherited from a d+ 1-dimensional gauge-fixing term used 
to facilitate the gauge-theoretic analogs of the manipulations of the previous 
sections. It turns out that the Feynman gauge is the simplest to use, and we 
will do so below. 

7 At least for so-called small gauge transformations, without topological properties; we 
will not discuss what happens with large gauge transformations in this paper. 




(27) 



A^UAfJJ-i + UdpU- 1 . 



(28) 
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The next question to discuss is that of periodicity of the gauge potential; 
in principle, periodicity need only be maintained up to a gauge transforma- 
tion. Because of the underlying gauge invariance of the basis used to calculate 



the trace in the partition function, one can (see, e.g., |L4| for thermal gauge 
theory and for the gauge-theory functional Schrodinger equation) intro- 
duce a projector which is an integral over all appropriate <i-dimensional gauge 
transformations into the functional integrals. Since small gauge transforma- 
tions are generated by 

fi{A(x)} = exp[i J d d xTr(KD i E, i )} (29) 

where E{ is the electric field (and i is a (^-dimensional index), one sees the 
well-known fact that gauge invariance is the same as imposing Gauss' law 
on the physical states occurring in the partition-function trace. As Ref. |H| 
shows, inserting the projector 

P = J[VA}tt{A(x)} (30) 

into the trace defining the partition function shows that the gauge partition 
function can be represented as a functional integral over all d+ 1 components 
of with strictly periodic boundary conditions, just as for the scalar case: 

A„(/3,x) = A M (0,x). (31) 

After inserting the projector P into the partition function, this quantity has 
a factor of the volume of the gauge transformations integrated over. The 
usual Faddeev-Popov gauge-fixing procedure must be applied to isolate this 
factor, subject to one proviso. These gauge transformations must obey a 
periodicity condition (which is simply periodicity of gauge- fixing ghosts): 

U(p,x) = U c U(0,x) (32) 

where JJq is an element of the center of the gauge group. Non-trivial elements 
of the center are important for discussing such non-perturbative phenomena 
as center vortices and confinement, but this will not be taken up here. 

Now we can show that using Feynman gauge allows for a trivial general- 
ization of the construction of the density matrix for scalars. In this gauge, 
the density-matrix propagator G v (t, x; r', x')^ is very simply related to the 
scalar propagator of equation (TO): 

G v (t, x; t', x% v = 5^ u G v (t, x; r', x'). (33) 
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Similarly, the gauge potential x) = which enters into the 

construction of the densitiy matrix analogously to the field <p(/3,x) of equa- 
tion (0), is just: 

A = fA^WLlMi (k) (34 ) 
A " J (2vr) d cosh(/3«; k /2) ^ W [6V 

where A^k) is the Fourier transform of the field A^x.) which is the argument 
of the density matrix. Repeating the arguments of Section [2.2| then shows 
that the thermal propagator for the gauge theory is again related to that of 
the scalar theory by a factor of 5^. 

The free term of the density matrix, constructed analogously to the scalar 
field density matrix (see equation (p0| )), is found to contain the Feynman 
gauge-fixing term in d dimensions, as expected. As mentioned above, it is 
also necessary to include periodic ghost contributions whose construction 
is entirely parallel to that of the scalar fields discussed above. Needless to 
say, the whole procedure can be carried out in any gauge, but certainly the 
simplest presentation of the technique is in the Feynman gauge. 

We will postpone applications of this formalism in gauge theories to future 
work, and now turn to a simple application in a scalar field theory. 



4 Applications to scalar theories 

The results of Section |2| provide us with a constructive method for deriving a 
DREA that can be interpreted as equivalent to the Landau "coarse-grained 
free energy" of the original microscopic theory. The T dependence of this 
free energy can be derived order by order. As long as we keep a cutoff (or use 
some other regularization method), we can proceed to compute correlations 
in the same manner as in any ordinary statistical mechanics problem. The 
renormalization required to obtain physical quantities can performed at the 
very last step, and will reflect the ultraviolet behavior of the original theory. 
On the other hand, the physical quantities themselves will be those of a d- 
dimensional theory whose statistical weights are dictated by the exponential 
of our reduced action. 

For practical applications, it is convenient to use the reduced theory not 
only because it reduces the number of dimensions, but also because it allows 
us to make contact with the vast literature on statistical mechanics. In this 
Section, we will illustrate how this comes about. 



14 



4.1 The DREA for ip 4 theory 

In order to compute correlations, we need not literally go through the steps 
described in Section 2 of adding sources J and integrating over both the r\ and 
(f) fields, followed by Legendre transformation, to actually find the DREA. A 
minor shortcut consists of obtaining the Feynman rules which emerge as we 
match the integrations with the rj integrations, resulting in a cancellation of 
ultraviolet divergences. For example, if we choose to compute (0(pi)0(p2)), 
pi = — p 2 = p, the combination of the rj and propagators which appear 
as internal lines in the various diagrams as a result of the integrations over 
those fields leads to a cancellation of divergences, and to the appearance of 
the usual thermal propagator G t h- As we will argue below, the resulting 
graphically-defined theory yields, for (0(pi)0(p2)), the sum of all graphs 
that would appear in a 4 theory in d- dimensions, with each internal line 
of momentum k corresponding to Gth(T2 — r i ; k) (the thermal propagator), 
and with each external line given by the G^(pj) propagator (j = 1,2) of 
equation ( pT|) multiplied by the factor cosh[w Pj (Tj — (3/2)]/ cosh(/3w Pj /2), all 
this integrated over T\ and r 2 . This is the field theory defined by the DREA. 

In order to show that the perturbative rules for the DREA are indeed 
the ones mentioned in the previous paragraph, we will consider the n-point 
correlation 

(0(pi) " " " 0(Pn)> = Z-\(3) J [V<j>] 0( Pl ) ■ ■ ■ 0(p n ) p[{3; 0, 0] , (35) 

where p\j3\ <fi, <fi] is given by ([□]). If we expand the interacting part, and make 
use of flTjD and (|H|) in the integral over rj, these rj-rj contractions will generate 
graphs with internal lines which correspond to G v , leaving the (p and <fi fields 
uncontracted. The final integral over <fr is to be performed with the Gaussian 
weight e Sp ^\ given by (^). This last integral will introduce (p-(p, and (p-<\> 
contractions. The former yield factors like those in ([22|), whereas the latter 
produce factors such as 

G t (v / 0Shla ^ ~,f " , (36) 
cosh(fjw Pj /2) 

for j = 1, . . . ,n. Now, the sum of (|22|) and (|23|) gives (|24l); to every rj-rj 
contraction there is a corresponding (p-(p contraction; and their weights are 
identical, as they come from the (p 4 graphs, with <p = (p + rj. All this shows 
that we shall have the graphs of 4 (x) theory, with G t ^ internal propagators, 
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and external lines given by with the r dependences of G t h and (35) 

integrated over. 

In the sequel, we will give an explicit example of how this works by exam- 
ining the case of broken (p A theory and performing a first order calculation. 



4.2 The phase transition in d = 2 

We will illustrate at one-loop level how the concept of dimensional reduction 
via the density matrix may be used to study the phase structure of ip A theory, 
as an example of a physical application. We will not be concerned with 
the fact that a one-loop calculation may not be an accurate description as 
we approach a phase transition, but will use this as an illustration of the 
techniques to compute physical masses in our scheme. Accurate or not, we 
will derive a criterion for a phase transition by demanding that the physical 
mass, the inverse of a correlation length, vanish at the transition point. 

The first nontrivial example we can treat is the case d = 2. The case d = 
is just quantum mechanics, its dimensional reduction leading to ordinary 
integrals, whereas the d = 1 case can be reduced to an effective quantum- 
mechanical problem, where no transition takes place. For d = 2 the theory 
is super-renormalizable, and only requires a mass renormalization. It has 
been studied in approximations involving both lattice computations [|T^. 
and analytic methods fl?| , I8|, [19 , pOf 
both. 



Our discussion will make contact with 



We start from an interacting Lagrangean defined by: 



(A/4!)(^-^-(l/2)mV 



(37) 



With real m the last term guarantees that the cutoff theory is in the broken 
symmetry phase. In order to derive the physical mass, we will compute the 
tadpole contribution to (0(p)0(— p)), and minimally subtract the equivalent 
term of the T = theory (ie, the tadpole in d+ 1 = 3). The self-energy, to 
first order, is given by: 



£ A (/3;p) = F(/3; p) 



A r d d k 



Gth(0;k) 



(38) 



F(flp) = 







2cosh^(/?w p /2) 



[1 + 



sinh(/3ty p 
(3w p 



(39) 
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where the last term of is the result of adding the tadpole with internal 
7] lines (which appears in the Appendix) to the one coming from contracting 
the first term in flP]) of the Appendix with the <fi propagator. In fact, the 
sum cancels the second term in (HTT) of the Appendix, restores r-translation 



invariance, and eliminates the d = 2 divergences explicitly shown in (fig, |49 



50) 



We now choose Ea(/5; p) = for (3 — > oo so that m is the physical mass 
parameter of the zero-temperature theory. Since F — > (l/w p ), we may infer 
the A-dependence to be attributed to \ipl in order to eliminate the T = 
divergence. Using the relation: 

' ' + ,.,,. t ,!„ TT » (40) 



2w k tanh(w k( 5/2) 2w k w k (e w ^ - 1) ' 



which splits the T = and T/0 parts of the last term of fl38|), and w^dw^ = 
kdk, we find that the divergence goes like A, whereas the finite remainder 
contributes: 

Sft = ^ln(l-e-/ r ). (41) 

The T-dependent correlation length £ = (m 2 + E^) -1 / 2 becomes infinite 
whenever T obeys: 

m 2 = -^ln(l-e- m / T ), (42) 
which, for T ^> m, becomes: 

AT T 

m 2 = ^ln(-)- (43) 
47r m 

There is always a critical temperature T c satisfying this equation for any 
positive value of A, with T c > m if A <C m. We emphasize that ( PB|) is 
exactly the same physical mass as occurs in the DREA as calculated from 
Z\fl; J], as claimed in [O . 

The discussion above has centered on rules for finding the DREA, or 
equivalently the matrix elements defined by the DREA. Once the DREA has 
been found, the whole theory is known and no further work is necessary. 
However, it is easier, in principle at least, to calculate the density matrix 
(or the corresponding action Sp of equation (|i~T|)) by only integrating over 
the r] fields, postponing the <ft integrations to a later step, such as numerical 
integration of the theory defined by Sd on a lattice. In such there will 
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remain spurious ultraviolet divergences, which must go away when matrix 
elements are computed via the functional integrals. This can be compared 
to a theory defined on a two-dimensional lattice. Its parameters depend on 
the physical T and on the cutoff (ie, the lattice spacing). We may then 
obtain the temperature where symmetry restoration occurs as a function of 
the cutoff. Removing the cutoff, so as to recover the continuum limit, will 
yield the physical transition T of the continuum theory. The calculation we 
have shown corresponds to a one-loop realization of this process; according 
to Refs. |L6|, 0, this is not terribly bad numerically. 

5 Conclusions 

We first introduced the notion of using the (negative logarithm of) the density 
matrix as defining an action Sd which would be a useful replacement for the 
usual DRFT, since it would not only be a theory in a reduced dimension 
(d + 1 — > d) but it would also be defined at all temperatures and might 
furnish a workable definition of the Landau- Ginzburg free energy, complete 
with information on phase transition temperatures coming from integrating 
over the Euclidean time dynamics. It then appeared that the action Sd could 
have spurious ultraviolet divergences not allowed in the partition function 
and its matrix elements, since these spurious divergences were absent in the 
T = theory. We showed that one could define a dimensionally-reduced 
effective action, the DREA, which was free of these spurious divergences as 
a result of cancellations between functional integrals over fields defining the 
Euclidean time evolution and functional integrals over the argument fields 
of Sd- The DREA is more difficult to calculate than the density matrix, 
since as an effective action it has infinitely-many terms involving all the 
Green's functions of the full theory. We have discussed here the graphical 
construction of the DREA correlations at the one-loop level. It is clear that 
dressed-loop techniques [|7|, [21]] apply for the DREA just as they do for any 
effective action; this allows for non-perturbative phenomena such as phase 
transitions to be investigated with the DREA. For example, corrections to 
the thermal propagator G t h(T2 — T i, k) can be formally summed so that the 
DREA is expressed entirely in terms of the dressed thermal propagator. 

We saw that one important difference between the DRFT and the DREA 
theory is that the DRFT may have ultraviolet divergences which are not in- 
herent in the finite-T version of the original d + 1 theory. These appear as 
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finite and calculable terms, generally involving In T, in the DREA theory and 
become the DRFT divergences at infinite temperature. We have studied an 
example of such divergences in the d = 2 + 1 <£> 4 model, studied by Einhorn 
and Jones |I8| using the DRFT. A one-loop mass divergence in the DRFT 
is replaced by a finite InT term in the DREA theory. A very similar phe- 
nomenon takes place in d — 2 + 1 non-Abelian gauge theory, whose DRFT is 
d = 2 gauge theory coupled to an adjoint scalar (the original time component 
of the gauge potential). This scalar in the DRFT has a logarithmic mass di- 
vergence which is actually a InT term at finite temperature, as calculated 
long ago by D'Hoker [E3]. The logarithmic divergence plays a very interesting 
role in understanding the transition from d = 2 + 1 gauge theory at zero T 
to the infinite-temperature limit, since the d = 2 + 1 theory has zero string 
tension for adjoint and similar Wilson loops but the corresponding DRFT 
has non- vanishing string tension for all representations of the Wilson loop. 



This question is under study by two of us |23[ . There will be other important 
uses of the DREA in investigations of d — 3 + 1 gauge theories as well (see, 
e.g., 0). 

We have emphasized that the DREA is capable of answering, in principle 
at least, questions about condensed- matter phase transitions which are not 
within the purview of the usual renormalizat ion-group approach to second- 
order phase transitions, such as the value of the phase transition temperature. 
Evidently, one aspect of the extension of the DREA to a dressed-loop expan- 
sion will continue to be the use of the renormalization group, so that one 
hopes to get both critical exponents and phase-transition temperatures from 
the DREA. Ultimately one hopes that the DREA can be truncated to rele- 
vant and marginal operators (for purposes of critical exponents) plus perhaps 
one or two more terms in the DREA expansion. 

Finally, we have noted but not explored at all the fact that the extension 
of the techniques used here for the density matrix can be straightforwardly 
extended to objects like the Wigner distribution function, which amounts 
to calculating off-diagonal matrix elements. This will be further explored in 
subsequent papers. 
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A Appendix 



Let us compute the contribution to the reduced action in first order for d = 2, 
in order to illustrate how divergences appear. From expression (|i"7|), we have: 



p x [(3\(j),(j)\ = --^po 



£ 4 (l) + 6 /G,(l; 1)^(1)+ 



(44) 



The last term in ( f44"D is independent of (p. The term in (p 2 corresponds to a 
tadpole graph which we will call T v . Then, using (|7|), (|8|) and ( |i~6|) we arrive 
at: 

_ r _fp_ 0(p)0(-p) f d*k_ 
Ir, ~J (2nY cosh 2 (/3 Wk /2) J {2-kY MP ' * j ' 

r , , s , , o, , nl N1 sinh(w k r) sinh[u> k (/? — r)l \ 

J„(p,k = / drcoshVp T-/3/2 1 ; . ryV^T ~" 46 

jo w k smh (p-u; k ) 

The integral over r is straightforward. It is convenient to use the identity: 

sinh(w k r) sinh[u; k (/3 - r)] = cosh 2 (/3w k /2) - cosh 2 [w k (r - /3/2)] , (47) 
and to split I v into a sum of three contributions, I' , I'' and J'". Explicitly: 

r = f3[l + smh((3w p )/((3w p )] 



4w k tanh(/3u> k /2) 



j// = 1 sinh[/?(w p + w k )] sinh[/?(w p - w k )] 

77 4-u; 2 . 8u' k (w p + u> k ) sinh(/?w k ) 8w k (w p — w k ) sinh(/3tt; k ) ' 

„, = g[l + sinh(/?Wp)/(/?w p )] 

" " 4w k sinh(/?u> k ) ' 1 ' 

hid = 2, we use w-^dw-^ = kdk. The integrals go from \m\ to (A 2 +m 2 ) 1 / 2 . The 
contribution from I' goes like A, the one from I" goes like In A, whereas the 
one from is finite as A — ► oo. In fact, the two contributions that diverge 
for large A correspond to the tadpoles in d = 3 and d = 2, respectively, if 
we use k 2 + m 2 as the inverse propagator. Although one might feel tempted 
to use these two types of subtractions to arrive at finite results, this is not 



what has to be done as shown in Section 2.2 
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